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(Abstract) This paper concerns with an economic analysis in e-learning environment. Regarding to the virtual nature of 
e-learning, uncertainties are inevitable. Learning consists of different integrated factors and due to real-time structures of 
e-learning in online content delivery, the components cost analysis help the decision makers to analyze the system economically. 
Here, the aim is proposing a comprehensive economic analysis model considering uncertainties in online course delivery. The 
input data being the online content delivery factors are considered to be fuzzy and the processing unit is a neural network based 
on Monte Carlo simulation as an adjustment core. The output of the proposed neural network is the cost estimation. The 
effectiveness and applicability of the proposed methodology is illustrated in a numerical example. 
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1. INTRODUCTION 

Despite the extensive development of the information 
technologies, information alone is not decisive to assure a 
competitive advantage. Global communication networks like 
the Internet enable the almost immediate exchange of 
information without limits of place and time [1]. However, it 
is not the information alone, but the knowledge and 
capability to assess, judge and apply these information in the 
process of a learning process [2]. Distance education is 
widely recognized as the alternative delivery system in 
which the student and the educator are separated either by 
distance or time or, in some cases, both. However, distance 
education and virtual learning are not new concepts, but 
rather have evolved from the ubiquitous correspondence 
courses of the past [3]. In the modern implementation, 
information or distributed learning technology is the likely 
connector between the learner, the instructor, and the 
offering site. Education at a distance may be as near as the 
on-campus residence halls or as far as a distant workplace 
[4]. 

Virtual learning is a term used to describe learning 
activities that can be accessed in a setting that is free from 
time and location constraints. Virtual learning comes in 
many forms, but it is about using computers to access 
information, courses and tools to develop skills and 
knowledge on all types of subjects. Virtual learning provides 
convenient and just-in-time access to learning programs and 
professional development materials. Virtual learning is 
learning that is enabled or supported by the use of digital 



tools and content [1]. It typically involves some form of 
interactivity, which may include online interaction between 
the learner and their teachers or peers. Virtual learning 
opportunities are usually accessed via the internet, though 
other technologies such as CD-ROM are also used in virtual 
learning. The last few years have seen a wave of digital tools 
and content being designed to facilitate the learning process 
[5]. Virtual learning enables different types of learning 
activities from those that rely on traditional teaching modes 
(e.g. the face-to-face lecture, tutorial or lab), and traditional 
media (e.g. books). It is increasingly used to support 
campus-based courses, and is bringing new dimensions to 
distance education [6]. 

Here, we propose a model to evaluate economic aspects of 
e-learning considering the uncertain parameters in online 
content delivery. Any e-learning is configured integrating a 
set of factors for online interaction. To evaluate the 
e-learning from cost view point, the factors should be 
assessed. Due to uncertainty in the status of the factors, an 
integrated system is designed to control the uncertainty and 
do the cost estimation computations. 

2. LITERATURE REVIEW 

For comparing the e-learning system and traditional learning 
system, two comprehensive studies were illustrated. 
Mahdavi et al. [7] compared traditional system with virtual 
educational system statistically in Iran. By the means of 
economical equations and statistical analysis they illustrated 
an in depth survey. They illustrated the best option for 



educational system is the combination of both systems. 
Fazlollahtabar and Sharma [8] compared traditional 
engineering educational system with the e-learning 
engineering educational system on the economic dimension 
using hypothesis testing approach in Iran. The comparison 
involved trend analysis and prediction based on costs and 
benefits of the two systems. Interestingly, the analysis 
revealed that the traditional system had greater advantage on 
the economic dimension. Several factors support the 
e-learning system despite the associated economic 
disadvantage. The final analysis provided results in favor of a 
blended system which takes advantage of the traditional and 
e-learning systems. 

Different studies have been worked out on cost 
optimization within e-learning environment. Mahdavi et al. [9] 
identified varied cost elements in e-learning educational 
system and optimized them by the means of mathematical 
programming. Then they proposed an effective method to 
estimate the learning cost between any two skills of learner 
using the grey relational analysis. Mahdavi et al. [10] 
developed their study combining the grey relational analysis 
and a radial basis function network to estimate the learning 
cost between any two skills after identification of varied cost 
elements in e-learning educational system and optimization 
by the means of mathematical programming. Fazlollahtabar 
and Yousefpoor [11] applied the cost elements in the 
e-learning educational systems and proposed a combination of 
grey relational analysis and a radial basis function network to 
estimate the learning cost between any two skills. An integer 
programming method was employed to demonstrate that it is 
possible to facilitate the acquisition of single skills by 
considering a set of useful compound skills. 
Fazlollahtabar and Mahdavi [12] proposed a neuro-fuzzy 
approach based on an evolutionary technique to obtain an 
optimal learning path for both instructor and learner. The 
neuro-fuzzy implementation helps to encode both structured 
and non-structured knowledge for the instructor. On the other 
hand, for learners, the neural network approach has been 
applied to make personalized curriculum profile based on 
individual learner requirements in a fuzzy environment. 

3. PROBLEM STATEMENT 

We consider an e-learning system for online content delivery 
to users. The components of such a system are hardware, 
software and human resource. Due to real-time structure of 
the system controlling decisions are accompanied with 
uncertainty. The major controlling decision in e-learning is 
economic view, cost estimation. Therefore, the costs and 
benefits of the system are treated using uncertain data being 
considered to be fuzzy. Since the estimation is performed 
using some sub-factors for hardware, software and human 
resource, a common and reliable method to relate sub-factors 
with main factors is neural network (NN). A configuration of 



neural network is shown in Figure 1 . 
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Figure 1 . A configuration of NN 

Also, we program the model for different time periods for the 
dynamism associated with the e-learning system. Thus, 
several outputs are obtained. To get unique cost estimation for 
the main components a Monte Carlo simulation is employed. 
The graphical presentation of the proposed methodology is 
shown in Figure 2. 



NN fa- 
estimation 



Hidden Layer 
input Layer Output Layer 



Fuzzy input 



I J 




Monte Carlo 
simulation for 
obtaining unique 
value 



Figure 2. The proposed methodology 



The basic materials are given next. 
4. BASIC MATERIALS 
4.1. Fuzzy Sets 

Here to imply the uncertainties in the VLEs, some basic 
definitions of fuzzy sets and fuzzy numbers are reviewed from 
[13, 14]. The basic definitions and notations will be used 
throughout the paper. 

Definition 1. A fuzzy set A in a universe of discourse X is 
characterized by a membership function /j^ ( f ) which 
associate©es with each element x in X a real number in the 
interval [0, 1 ] . The function value ju^ {x ) is termed the grade 
of membership of x in A (Buckley, 1987). 
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Definition 2. A fuzzy number A is a fuzzy convex subset of 
the real line satisfying the following conditions: 

(a) jU^ (X ] is piecewise continuous; 

(b) JU^ [X j is normalized, that is, there exists fflG?i with 
jil^ {m. ) = 1 , where m is called the mean value of A [15] . 
Definition 3. A triangular fuzzy number a can be defined by 
a triplet ( ) . Its conceptual schema and 
mathematical form are shown by Equation (1): 



0, 


X <dy 


x -a { 


a { <x <a 2 




a 2 -a j 




a 3 -x 


a 2 < x <a 3 


a 3 -a 2 ' 


0, 


a } < x . 



(1) 



A triangular fuzzy number a in the universe of discourse X 
that conforms to this definition has been shown in Figure 3. 



Ms ( x ) 




a 1 a 2 a 3 X 

Figure 3. A triangular fuzzy number a . 
Definition 4. Assuming that both d — (a l ,a 2 ,a 3 ) and 

b=(b l ,b 2 ,b 3 ) are triangular numbers, then the basic 
fuzzy operations are: 



dxb =(a l y.b x ,a 2 xb 2 ,a 3 xb 3 ) 



for multiplication (2) 



a +b =(a x +b 1 ,a 2 +b 2 ,a 3 +b 3 ) for addition (3) 



4.2. The Backpropagation Neural Network 



The backpropagation algorithm trains a given feed-forward 
multilayer neural network for a given set of input patterns 
with known classifications. When each entry of the sample set 
is presented to the network, the network examines its output 
response to the sample input pattern. The output response is 
then compared to the known and desired output and the error 
value is calculated. Based on the error, the connection weights 
are adjusted. The backpropagation algorithm is based on 
Widrow-Hoff delta learning rule in which the weight 
adjustment is done through mean square error of the output 
response to the sample input. 
Algorithm 1: General steps of backpropagation. 

1. Propagate inputs forward in the usual way, i.e., all outputs 
are computed using sigmoid thresholding of the inner product 
of the corresponding weight and input vectors. All outputs at 
stage n are connected to all the inputs at stage n+l 

2. Propagate the errors backwards by apportioning them to 
each unit according to the amount of the error the unit is 
responsible for. 

We now discuss how to develop the stochastic 
backpropagation algorithm for the general case. The 
derivation is simple, but unfortunately the book-keeping is a 
little messy. The following notations and definitions are 
needed: 

Xj : input vector for unit j (x ; , = ith input to the y'th unit) 
Wj : weight vector for unit j (wp = weight on xj,) 
Zj = Wj .Xj : the weighted sum of inputs for unit j 
Oj : output of unit j ( O j = <j(Zj ) ) 
tj : target for unit j 

Downstreamij) : set of units whose immediate inputs include 
the output of j 

Output : Set of output units in the final layer. 
Since we update after each training example, we can simplify 
the notation somewhat by assuming that the training set 
consists of exactly one example and so the error can simply 
be denoted by E. 

8E 



We want to calculate 



corresponding to each input 



Definition 5. The D p q -distance, indexed by parameters 
1 < p < co and < q < 1 t between two fuzzy numbers 

CI and b is a nonnegative function given by: 



(l-q)^ja a -b-\ P da + q^jal-b^ da " , p<x>, (4) 
(l-q) sup (I a~ -b~\ ) + q inf (la* -6*1 ] , p = oo. 



weight Wji of each output unit j. Note first that since Zj is a 
function of Wp regardless of where in the network unit j is 
located, 



dE _ dE dzj _ dE 



dw ji d Zj ^.a d Zj 



8E 



.x. 



(5) 



Furthermore, zfL is the same regardless of which input 

8 Zj 

weight of unit j we are trying to update. So, we denote this 
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quantity by 8 ■ . 

Consider the case when j is an output unit. We know that 
E = \ Z(h-v(zJ) 2 - ( 6 ) 

Since the outputs of all units A: ^ j are independent of w 7 „ 
we can then drop the summation and consider just the 
contribution to E by j and we call it 8 ■ : 



SE 5 1, , 2 , , do j , x 8 , , rT , 

8, =3-=^ ~°> f - (t > -^t^' -°' ) ^ a(z ' ) (?) 



we obtain: 

S k =Oj<X-Oj) 



= -(tj -o J )(l-a(z J ))<7(z J ) = -(tj -o j )(\-o j )o j . 



Thus, 



dE 



Aw ,, = -77 = rj8,x , . 



(8) 



Now, consider the case when j is a hidden unit. Like before, 
we make the following two important observations: 

1 . For each unit k downstream from j, Zk is a function of Zj. 

2. The contribution to error by all units I ^ j ' , in the same 
layer as j, is independent of w jt . 

8E 

We want to calculate for each input weight w /; for each 

dw jt 

hidden unit j. Note that w,, influences just Zj which influences 
oj which influences Zh Vk e Downstream(j), each of which 
influences E. So, we can write, 

dE £ afi&^d ^^, ^ gg 8z t 8oj x -(9) 

SWji k„DownstreamU) ^0 j 8Zj 8w jt teD„wnst m ,m(J) 8z t do j dZj 

Again, note that all the terms except x ; , in (9) are the same 
regardless of which input weight of unit j we are trying to 

update. Like before, we denote this common quantity by 8 . 



dE c dz k do; 
Also, note that 8 t , W, ; and = 0,(1-0,). 

dz k doj k] dij J 1 

Substituting them in (7), 

£ dE_8^do L= ^ s (10) 

keDownstream(j) j j k&Downstream{ j) 



£4 

keDownstream ( 7 ) 



(11) 



4.3. Fuzzy neural networks 

Fuzzy neural networks (FNN) have been proposed as a 
knowledge engineering technique and used for various 
application domains by several authors including Yamakawa 
et al. [16]. Fuzzy sets were initially introduced into neural 
networks in 1993 [17], but has not experienced much 
development until recently. Since 1995, much attention has 
been focused on incorporating fuzzy reasoning into neural 
networks [18]. A neural network is considered a fuzzy neural 
network if the signals and/or the weights in the system are 
based around fuzzy sets [19]. The fuzzy neural network 
systems (Neuro-fuzzy systems) combine the advantages of 
fuzzy logic systems and neural networks have become a very 
active subject in many scientific and engineering areas, such 
as, model reference control problems, PID controller tuning, 
signal processing, etc. The fuzzy neural network (FNN) 
system is one kind of fuzzy inference system in neural 
network structure [20, 21, 22, 23]. 

4.4. Monte Carlo simulation 

Monte Carlo simulation is a comprehensive approach for 
analyzing the behavior of some activities, plans or processes 
that involve uncertainty. If we face uncertain or variable 
market demand, fluctuating costs, variation in a manufacturing 
process, or effects of weather on operations, or stochastic 
activity time we can benefit from using Monte Carlo 
simulation to understand the impact of uncertainty, and 
develop plans to mitigate or otherwise cope with risk. 
Whenever we need to make an estimate, forecast or decision 
where there is significant uncertainty, we'd be well advised to 
consider Monte Carlo simulation [24]. 

Monte Carlo simulation is a method for iteratively 
evaluating a deterministic model using sets of random 
numbers as inputs. This method is often used when the model 
is complex, nonlinear, or involves more than just a couple 
uncertain parameters. The Monte Carlo method is just one of 
many methods for analyzing uncertainty propagation, where 
the goal is to determine how random variation, lack of 
knowledge, or error affects the sensitivity, performance, or 
reliability of the system that is being modeled. Monte Carlo 
simulation is categorized as a sampling method because the 
inputs are randomly generated from probability distributions 
to simulate the process of sampling from an actual population. 
So, we try to choose a distribution for the inputs that most 
closely matches data we already have, or best represents our 
current state of knowledge. The data generated from the 
simulation can be represented as probability distributions (or 
histograms) or converted to error bars, reliability predictions, 
tolerance zones, and confidence intervals. 

The following five steps are proposed to implement Monte 



Carlo, 

Step 1: Create a parametric model, y = f{x h x 2 , x q ). 

Step 2: Generate a set of random inputs, x u , x i2 , x iq . 

Step 3: Evaluate the model and store the results as y,. 

Step 4: Repeat steps 2 and 3 for = 1 to n. 

Step 5: Analyze the results using histograms, summary 

statistics, confidence intervals, etc. 

5. PROBLEM FORMULATION 

Here, we formulate our proposed model including uncertainty 

in economic evaluation of e-learning online content delivery 

factors. Initially, we determine the effective factors in online 

content delivery to be hardware, software and human resource. 

Learning in e-learning is an integrated network of some 

factors. Learning depends on the suitable functioning of 

hardware (servers, computers, ...), software (learning 

management systems), and human (instructor, 

administrator, ...). Each of these factors could be assessed 

economically. While the status of these factors associated with 

uncertainty, we need a method to evaluate them during time. 

To identify whether a factor is in its appropriate status or not, 

two economic expressions are applied. Net present value and 

benefit to cost ratio are two economic expressions that are 

used to evaluate the economic status of a factor. 

The economic model is as follows: 

r. rate of return 

/: unit of time 

B,: benefit at time t 

C t : cost at time t 

NPV: net present value 

BCR: benefit to cost ratio 

To gain the net present value we use, 

(B,-C t ) 



npv = y 

m (! + >-)' 
In fuzzy case we have, 

NPV = Y V — r 1 

To obtain benefit to cost ratio we use, 

f to) 

^ CM' 

In fuzzy case we have, 



(12) 



(13) 



(14) 



M 
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(15) 



f gj 

These equations are used to estimate the cost and benefit in a 
neural network proposed in Section 4.2. To do the 
computations in fuzzy environment for all of the factors we 
use equations (13) and (15). Consider that all costs, benefits, 
and rate of returns are uncertain and triangular fuzzy numbers. 
As we described in Section 4. 1, we do fuzzy basic operations 

to determine NPV and BCR for each factor. After that, 

we have to compare the resulted NPV and BCR to gain 

the optimal (maximum) one. This means that it is necessary to 

have a method for ranking or comparing fuzzy numbers. An 

operator ■< for ordering fuzzy numbers can be defined as 

follows: 

A<B (a 1 <fc 1 )A(a 2 <fc 2 )A(a 3 <fc 3 )A(a 4 <fc 4 ). (16) 
However, this relation is not a complete ordering, as fuzzy 



numbers A , B satisfying 
3 i, j e {1, 2, 3, 4}) : (a,. <b t )A (a,. > b, ) 



(17) 



are not comparable by ■< . 

Here, we introduce a new fuzzy ranking method for fuzzy 
numbers. Let us consider fuzzy max operation defined in the 
following way: 



Max value(5,&) = (max(a 1 ,ft 1 ),max(fl 2 ,fe 2 ), 
max(a 3 ,& 3 ), max(a 4 ,& 4 )). 



(18) 



It is evident that, for non-comparable fuzzy numbers a, b , the 

fuzzy max operation results in a fuzzy number different from 
both of them. To alleviate this drawback, we propose a new 
method based on the distance between fuzzy numbers. We use 
the distance function introduced in [25]. The analytical 

properties of D ] p q depend on the first parameter p, while the 

second parameter q of D p q characterizes the subjective 

weight attributed to the end points of the support; i.e., ( a + a , a~ ) 
of the fuzzy numbers. If there is no reason for distinguishing 



any 



side of the fuzzy numbers, D l j s recommended. 



"'2 



Having q close to 1 results in considering the right side of the 
support of the fuzzy numbers more favorably. Since the 
significance of the end points of the support of the fuzzy 

_ 1 

numbers is assumed to be the same, then we consider — — . 
For triangular fuzzy numbers 

a=(a l , a 2 , fl 3 ) andb = , Z? 2 , ^3 ) , the above distance 
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To illustrate the effectiveness and applicability of the 
proposed methodology a short numerical example is 
illustrated. Input data are collected from the experts for the 
components and inserted to the NN for cost estimation 
purpose. A feedforward backpropagation network can be used 
to approximate a function to relate the sub-factors to the main 
factors. To facilitate the computations of backpropagation 
neural network, MATLAB 7 user interface, NNtool, was 
applied. A feedforward network was programmed with 
sub-factor inputs, ten hidden units with logistics activation 
function, and the outputs. A configuration of the feedforward 
neural network is shown in Figure 4. 



with p — 2 and Q — — is then calculated as: 



D.a.b 



i=l fe{l,2) 



(19) 



* IW{1,1} 

bin 




L 



LW{2,1} 



b{2} 




<3t 



10 



Thus, we can acquire the distance between CL,b and 

M a V using the proposed distance function. We proposed the 

following algorithm for the comparison of fuzzy numbers. 
Algorithm 2: Comparing two triangular fuzzy numbers 

Input: Two fuzzy numbers ( L- , i = 1, 2 ) 

Output: Maximum between two triangular fuzzy numbers 

(L max ) 

Step 1: Calculate the maximum value (M V ) by Equation 
(8); ( M a V stands for Max Value). 

Step 2: Find the distance of M a V from (L ( , i =1, 2 ) 
using Equation (19). 

Step 3: Determine L mdx with the highest distance. 

Therefore, we gain several estimations for component and 
using Monte Carlo simulation a unique value is in hand. 



Figure 4. A configuration of the feedforward neural network 

Using NNtool, the data were inserted and the required settings 
were made to train the data to obtain an appropriate pattern. 
An applied user interface is presented in Figure 5. 



Network: networkl 
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rnujnc 
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show 

time 



%1 Train Network 



Figure 5. A user interface in MATLAB 

Note that, the input data were triangular fuzzy numbers to 
imply the uncertainty of the data. Now, we apply Monte Carlo 
simulation to obtain the unique value for each main 
component as explained before. According to a sampling 
procedure of occurrences, we compute the variable cost 
corresponding the cumulative probability of occurrence for 
each component which is reported in Table 1 . 



6. NUMERICAL EXAMPLE 
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Table 1. The variable times corresponding the probability of occurrence 



Hardware 


Software 


Human resource 


Cumulative Probability 


Sorted times 


Cumulative Probability 


Sorted times 


Cumulative Probability 


Sorted times 


0.20 


3.8 


0.15 


9.6 


0.08 


1.65 


0.45 


3.9 


0.28 


9.7 


0.026 


1.85 


0.75 


4 


0.53 


9.9 


0.48 


2 


0.85 


4.1 


0.73 


10 


0.72 


2.1 


0.93 


4.2 


0.9 


10.3 


0.88 


2.2 


1 


4.25 


1 


10.35 


1 


2.25 















Here, according to the Monte Carlo process we generate five 
series of random numbers as given in Table 2. 

Table 2. The generated random numbers 





HARDWARE 


SOFTWARE 


HUMAN RESOURCE 


1 


0.632175 


0.642 


0.381758 


2 


0.238231 


0.0211 


0.446354 


3 


0.567427 


0.1123 


0.278362 


4 


0.613292 


.78532 


0.872533 


5 


0.013479 


0.34578 


0.253833 



Now, we accommodate the random numbers of each 
component with the corresponding component estimation in 
Table 1. The following component cost estimation, Table 3, is 
obtained for the components associated with the average cost 
for each component of e-learning. 

Table 3. The obtained activity times 





HARDWARE 


SOFTWARE 


HUMAN 
RESOURCE 


1 


4 


10 


2 


2 


3.9 


9.6 


2 


3 


4 


9.6 


2 


4 


4 


10.3 


2.2 


5 


3.8 


9.9 


1.85 


Average 


3.94 







Considering the aggregation rule, we have the following cost 
estimation: 

Total estimated cost= 3.94+9.88+2.01= 15.83, 

Then, the total estimated cost equals to 15.83 unit of money. 

7. CONCLUSIONS 



Here, we proposed an integrated methodology based on fuzzy, 
neural network and Monte Carlo simulation to analyze the 
economic status of uncertain factors of e-learning. For 
economic assessment net present value and benefit to cost 
ratio were applied. All factors were considered triangular 
fuzzy numbers and using economic equations and fuzzy basic 
operations the neural network performed the cost estimation 
and fuzzy values as the output were obtained in different time 
periods. To compare those fuzzy numbers an algorithm was 
introduced and the ranking of e-learning factors is gained. The 
Monte Carlo simulation was employed to aggregate the time 
based dynamic values and obtained unique cost estimation for 
each component for online content delivery in e-learning. The 
applicability of the methodology was illustrated in a 
numerical study. 
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